On the inertia of heat 
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Does heat have inertia? This question is at the core of a long-standing controversy on Eckart's 
dissipative relativistic hydrodynamics. Here I show that the troublesome inertial term in Eckart's 
heat flux arises only if one insists on defining thermal diffusivity as a spacetime constant. I argue 
that this is not the most natural definition, and that all confusion disappears if one considers instead 
the space-dependent comovmg diffusivity, in line with the fact that, in the presence of gravity, space 
is an inhomogeneous medium. 



Eckart was the first to consider the problem of heat 
conduction in general relativity. In the seminal 1940 pa- 
per he proposed as that the general-relativistic the 
heat equation 



u a V a T + 9T + V a q a = 0. 



(1) 



Here u a is the 4-velocity of the thermal medium, 9 = 
V a u a its expansion, and q a the heat flux0 To close this 
equation, he proposed the constitutive relation 



q b = -n (D b T + a b T) 



(2) 



where D b T = (u a u b + 5 c b l )V a T is the spatial gradient of 
T (in the local rest frame associated to u"), a b = u a V a u b 
the local acceleration, and kq the thermal diffusivity. 
This ansatz, he argued, is the simplest one that is consis- 
tent with the second law of hydrodynamics. 

In spite of early criticisms of Eckart's approach 0, 
the system ([I])-® (or a variant thereof due to Landau and 
Lifshitz was considered the standard formulation of 
relativistic thermal dynamics until the late 70's. Together 
with the relativistic Navier-Stokes equations, it formed 
the theory of "dissipative relativistic hydrodynamics" ; its 
connection to relativistic kinetic theory was believed to 
be well understood |||. 

The situation changed drastically with two major con- 
tributions in the field: Israel and Stewart's development 
of a "second-order" formalism bridging between the ki- 
netic (microscopic) and hydrodynamic (macroscopic) lev- 
els of description of relativistic transport phenomena 
[H, 0]; and the discovery by Hiscock and Lindblom of 
"generic instabilities" in Eckart's "first-order" formalism 
(In this context, the expressions "first order" and 



"second order" refer to the parabolic and hyperbolic na- 
ture of the corresponding system of partial differential 
equations.) The consensus soon became that Eckart's 
theory is both acausal and unstable, and thus "unaccept- 
able". Given the high complexity of the Israel-Stewart 
"extended thermodynamics" , with involves many more 
dynamical fields than the standard hydrodynamic ones 
(temperature, velocity, pressure), the field has remained 
— to this day! — without an accepted theory of relativis- 
tic dissipative phenomena. This is all the more frustrat- 
ing as new, important applications have arisen in high- 
energy physics, in particular regarding the dynamics of 
the quark-gluon plasma. 

The bottom line of the recent criticisms to Eckart's 
theory is the presence of a temperature-acceleration cou- 
pling a b T in the constitutive relation. This term, coined 
the "inertia of heat" , has no analogue in Fourier's con- 
stitutive relation, and is allegedly responsible for the un- 
physical behavior of the heat equation Among the 
old and recent arguments to this effect, one reads that: 
it is not consistent with the relativistic version of the 
first principle it depends on an arbitrary choice of 
frame (lo| : it is not a gradient, while linear irreversible 
thermodynamics requires that all dissipative fluxes are 
gradients [111 ]: heat is a non-mechanical form of energy, 



so it should not have "inertia" I12H14H ; it is not consistent 
with relativistic kinetic theory [ll| ; it is the source of the 
Hiscock-Lindblom "generic instabilities" fl5, 16|. 

My purpose in this note is not to address these ob- 
jections head-onH Rather, I wish to discuss the physical 
origin of the inertia of heat. But before that, let me make 
three preliminary observations. 

First, the issue of causality is not an issue at all. The 
usual Fourier heat equation also propagates signals faster 



1 Equivalcntly, this equation can be written as the continuity equa- 
tion V a p a = for the energy flux p a = Tu a + q a (in units where 
the heat capacity is 1). 



2 The issue of the consistency of Eckart's ansatz with relativistic 
kinetic theory will be addressed elsewhere |17|| . 



2 



than light, in patent contradiction with the much slower 
values of microscopic thermal velocities. All this means is 
that we should not forget that dissipativc hydrodynamics 
is a long-time, large-scale approximation of kinetic the- 
ory, which breaks down when spacetime gradients arise 
on scales comparable to the mean free flight and mean 
free path. In its regime of validity, the heat equation is 
perfectly consistent with causality, and the same holds for 
Eckart's equation. This point was made repeatedly in the 
literature pUGl], 

and should by now be widely acknowl- 
edged. Furthermore, it has been showed that alternative 
"hyperbolic theories" of the Israel-Stewart type, which in- 
volve more dynamical variables and free functions, gener- 
ally have the same physical content as parabolic theories 
a la Eckart 0. 

Second, the "generic instabilities" discovered by His- 
cock and Lindblom @, Q and referred to in 15, lg per- 
tain to the coupled heat-fluid dynamics, i.e. the relativis- 
tic Fourier- Navier-Stokes system. The issue at stake here 
is whether it is consistent to include Eckart's heat flux 
(J2J) into the stress-energy-momentum tensor T a b entering 
the hydrodynamical equation V a T ab = for a dissipative 
fluid; some authors have argued that it should not [l2j ]. 
Here, I would like to focus instead on the heat dynam- 
ics itself, assuming a prescribed background fluid velocity 
u a . This is as sensible approximation as Fourier's theory 
for heat conduction, and is valid in a regime where the 
backrcaction of the heat flux on the medium's dynamics 
is negligible. In this approximation, the no-go result of 
[1, E|, however important in itself, is not relevant; in par- 
ticular it must not prevent us from considering the issue 
of the inertia of heat in itself. 

Third, there is an easy way to understand the meaning 
of the inertia of heat: by considering thermal equilibrium 
in static spacetimes. Denoting £ a a timelike Killing vec- 
tor and x = ( — £ a £a,y/ 2 the redshift factor, and observ- 
ing that the acceleration of the hydrostatic congruence 
u a _ £<y^ i s ab = Dbx/Xi we see that gt, = is equiva- 
lent to xT = const, as prescribed by Tolman's law [20[ . In 
other words, the inertial term in Eckart's ansatz is exactly 
what it takes to ensure that the equilibrium temperature 
distribution compensates gravitational redshifts. This ar- 
gument is very compelling; in fact, it was almost always 
cited in the early reviews of relativistic hydrodynamics, 
as e.g. in 

I now come to the thesis of this note: the disturbing 
inertial term in Eckart 's ansatz comes from an awkward 
interpretation of the notion of "thermal diffusivity". Let 
me explain. 

In Eckart's approach, and in all subsequent work I am 
aware of, thermal diffusivity is defined as a constant co- 



efficient Ko, as in ©. However, dimensional analysis 
teaches us that diffusivity has dimensions (length) 2 /time, 
and this implies that its value should change if the unit 
of time changes. Now, in a curved spacetime, this is pre- 
cisely what happens from point to point. 

For clarity, assume as before that the thermal medium 
is in hydrostatic equilibrium, so that its instantaneous 
rest frame can be described by static comoving coordi- 
nates (t,x l ) such thai|f| 



ds 2 = —x{ x ) 2 dt 2 + qijdx t dx-' . 



(3) 



Here x is again the redshift factor and qij is the comoving 
spatial metric, with associated covariant derivative Di. In 
these coordinates, Eckart's heat equation become 

X - 1 d t T = Ko (D i + a i )(D i T + a i ), (4) 



d t T — k x(AT + 2a i D i T + a 2 T), (5) 

where A = DiD 1 is the spatial Laplacian. Now, defining 
the function 



k{x) = Kqx(x) 
it is easy to see that ([5]) is equivalent to 
d t T = A(kT). 



(6) 



(7) 



This is an interesting observation: Eckart's equation 
neatly simplifies when we use the function k(x) instead 
of the constant kq. In particular, all terms involving the 
acceleration a 1 disappear. But what is the interpretation 
of k{x)1 

My answer is that n(x) should be interpreted as the 
comoving diffusivity, and heat conduction in general rel- 
ativity as diffusion in a inhomogeneous medium. 

The most useful definition of diffusivity is in terms 
of the mean square displacement of a Brownian parti- 
cle whose probability density satisfies the corresponding 
heat equation. If (d 2 (x s ,xo)) is the mean of the squared 
spatial distance between the position x s of a Brownian 
particle at time s and its starting point xo, we can de- 
fine the diffusivity as lim. s ->o(d (x s , xq))/6s. Now, from 
(JT)), we find that if s is the proper time r along a static 
worldlinc initiated at Xq, then 



lim 



{d 2 (x T ,x )} 
6r 



«0) 



(S) 



3 We assume that u a is irrotational. 
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while if s is the comoving time coordinate t, we have 



lim 



(d 2 (x t ,x )) 
6t 



= k(Xq). 



(9) 



This observation justifies the names proper diffusivity for 
kq and comoving diffusivity for k{x). The point is that, 
while the former is constant in space, the latter is not 
(unless there is no gravity). This property is precisely 
what makes the comoving diffusivity physically relevant: 
it expresses the fact that, if time runs faster in some re- 
gion of space, then diffusion is also faster in these regions. 
This is the origin of the "inertia of heat" . 

Of course, the notion of proper diffusivity is also useful, 
since it is the one related to the proper time measured by 
physical clocks. The point is that this invariant notion 
is at variance with an intuition based on the notion of 
space (as opposed to spacetime). Indeed, this notion of 
involves a foliation of spacetime indexed by the comoving 
time coordinate t, and therefore requires that diffusivity 
be measured in units of this global time - where it is not 
homogeneous but space-dependent. Thus, I propose that 
the unease with Eckart's inertial term comes from a clash 
between spatial intuitions and the covariant formulation 
of the heat equation JT]) . 

Let me now step back and reconsider the derivation of 
the general-rclativistic heat equation from this perspec- 
tive. In comoving coordinates, the thermal medium is 
inhomogeneous, with a space-dependent diffusivity n(x). 
As usual, we may write the continuity equation in the 
form 



dtT = -DiQ 1 



(10) 



with Q l the comoving heat flux, involving the gradient of 
k and T. Now, given that diffusivity is now a function 
of space, the form that Q % should take is not obvious at 
priori: should it take the "Fourier" form 



U = -kD*T 
or the "Fokker-Planck" form 

Q* = -D\kT), 



(11) 



(12) 



"microscopic bias" implying that the mean position (x s ) 
of Brownian particles changes with time, the correct form 
is the Fokker-Planck one |22| • If we now remember that 
heat is carried by the very same Brownian particles which 
make up the thermal medium, and that in the comoving 



frame their mean velocity vanishes, we see that the right 
pick for the heat equation is (|12p . This gives ([7]), which 
we saw is indeed Eckart's equation (in hydrostatic equi- 
librium). Observe that (|12[) is a gradient, as required by 
linear irreversible thermodynamics. 

With hindsight, it is not particularly surprising that 
the relativistic heat equation is more easily interpreted in 
the comoving frame than in covariant form ([T]), in spite of 
the fact that the corresponding time coordinate t is not 
observable: dissipative phenomena are by nature alien to 
covariance. They are associated to the production of en- 
tropy, viz. with a thcrmodynamical arrow of time. Now, 
it so happens that the natural time variable associated to 
this arrow is not the local proper time r: it is the global 
comoving time t. 
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